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ABSTRACT 

The statistics of peaks in weak gravitational leasing maps is a promising technique to 
constrain cosmological parameters in present and future surveys. Here we investigate 
its power when using general extreme value statistics which is very sensitive to the 
exponential tail of the halo mass function. 

To this end, we use an analytic method to quantify the number of weak leasing peaks 
caused by galaxy clusters, large-scale structures and observational noise. Doing so, we 
further improve the method in the regime of high signal-to-noise ratios dominated by 
non-linear structures by accounting for the embedding of those counts into the sur¬ 
rounding shear caused by large scale structures. 

We derive the extreme value and order statistics for both over-densities (positive 
peaks) and under-densities (negative peaks) and provide an optimized criterion to 
split a wide held survey into sub-helds in order to sample the distribution of extreme 
values such that the expected objects causing the largest signals are mostly due to 
galaxy clusters. We hnd good agreement of our model predictions with a ray-tracing 
A-body simulation. 

For a Euclid-like survey, we hnd tight constraints on erg and Dm with relative uncer¬ 
tainties of ~ 10“^. In contrast, the equation of state parameter wq can be constrained 
only with a 10% level, and Wg, is out of reach even if we include redshift information. 

Key words: cosmology: theory - gravitational leasing: weak - cosmological parame¬ 
ters - large-scale structure of the Universe 


1 INTRODUCTION 


The cluster mass function has proven to be a powerful tool to 
constrain cosmological parameters (e.g. Haiman et al. l200T| 
Majumdar fc Mohr||2004| |Wang et al.| 2004| and |Vikhliniii| 


et al. 


20091. However, the mass is not a direct observable 


and can only be inferred indirectly. Not even gravitational 
lensing provides a direct measure of mass because it only 
constrains the gravitational potential, which is a local quan¬ 
tity in contrast to the mass. It is therefore advantageous to 
replace the mass with a well defined observable quantity to 
constrain cosmological parameters in an unambiguous way. 
Additionally, this observable needs to be predictable from 
theory in a reliable and clear way. There are different ways 
to tackle this problem. For instance the mass can be replaced 


with the X-ray temperature of galaxy clusters (Angrick & 
Bartelmann|200^ Angrick et al.|2015 l. 


In this paper we replace the mass with the signal-to- 
noise {S/N henceforth) ratio of weak lensing shear detec¬ 
tions ( Maturi et al.||2009[ 201 1| which is based on the work 
of Bardeen et al.||lS)^ and van Waerbeke|2()00 \. We then ap¬ 
ply general extreme value (GEV) and order statistics to the 
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abundance of S/A-ratio peaks in cosmic shear maps to re¬ 
produce the distribution of the highest detections (i.e. with 
the highest S/N) correctly. Cosmic shear peaks originate 
from three sources: i) From noise due to the intrinsic ellip- 
ticity distribution and the finite number of sources, ii) from 
chance alignments of large-scale structures (LSS) and iii) 
from cluster-sized dark matter halos. The highest peaks in a 
weak lensing map will be due to clusters. Therefore they will 
carry similar cosmological information as the mass function. 
Given the constraining power of weak lensing peak counts 
( Marian et al.|[2C)09 Kratochvil et al. 2010 and Dietrich & 
Hartlap 20101, the abundance of weak lensing detections has 
been predicted and used in|Maturi et al.|||2009p , |Maturi et al. 
(20111, Cardone et al. (20131, Cardone et al. (20141, Petri 


et al. (20141, Mainini & Romano (20141 and Lin & Kilbinger 


Under-densities (voids) produce weak lensing signals as 
well as over-densities. It is possible to model the abundance 
of negative detections by a Gaussian random field, as those 
objects evolve mostly linearly. In contrast the over-densities 
produce detections which are due to halos and therefore due 
to highly non-linear objects. 

We use the abundance of detections of both positive 
(over-densities) and negative (voids) to extract their respec- 
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tive GEV and order distributions. This was already done for 


the mass function in Davis et al. (20111, Waizmann et al. 
(20111, Waizmann et al. (20121 and Waizmann et al. 120131. 


This approach has the advantage that only the highest weak 
leasing peaks are taken into account which are prominent 
features in weak lensing maps. Even though information is 
in principle lost when considering the highest peaks only, 
these objects will populate the exponential tail of the mass 
function and are therefore very sensitive to variations in the 
cosmological parameters. Furthermore this approach avoids 
the assumption of a Gaussian likelihood and links the ob¬ 
servable directly to its distribution. 

We focus on the statistical distribution of the highest 
positive and negative weak lensing peaks, which trace the 
most massive objects and the most under-dense regions re¬ 
spectively. This approach has some advantages: First iden¬ 
tifying weak lensing detections, i.e. connected regions that 
have a weak lensing signal larger than some threshold, al¬ 
lows to avoid the poorly constrained mass of a cluster as 
an observable, (cf. Maturi et al. 2009 20111. Furthermore 
the physical assumptions entering the model are minimized 
since the theoretical prediction is directly linked to observa¬ 
tional data. Counting detections above different thresholds 
allows to construct an observable similar to the mass func¬ 
tion. However, the S/Wthresholds define detections in the 
hrst place. Therefore measurement errors are mainly due to 
hnite number of detections. Second, by focusing on the high¬ 
est peaks only the most prominent weak lensing peaks are 
important, which are defined very well. This implies that 
the observable is very well defined and its distribution is 
also known directly from the model. Thus any assumption 
about the distribution can be omitted 

We investigate the impact of changes in the cosmologi¬ 
cal parameters and compare the resulting distributions are 
compared to numerical A'^-body simulations carried out by 


Borgani et al. (2004) using Gadget-2 (Springel (2005). We 


then investigate the constraining power (with respect to the 
matter density flm, the amplitude of the matter power spec¬ 
trum (T8 and the equation of state parameter wo for spatially 
flat cosmologies) of this method by performing a maximum- 
likelihood analysis on mock data with the characteristics 
expected for future space-based sur veys, such as for exam - 
ples the Euclic0ESA space mission Laureijs et al. (2011), 
for which we assume a sky coverage of 15.000 deg^ and a 
galaxy number density of 40arcmin“^. 

Throughout this work we use the cosmological base pa¬ 
rameters from |Planck-Collaboration et al.| ( |2013[ ) , namely 
the matter density parameter = 0.314, the cosmologi¬ 
cal constant density parameter Ha = 1 — the baryon 
density parameter Q.\,oh? — 0.02206, the Hubble constant 
h = 0.686 and the normalization of the matter power spec¬ 
trum (Tg = 0.834 as hducial values. 

The structure of the paper is as follows: in Sect. we 
summarise the theoretical model used for the abundance of 
cosmic shear peaks describing the different contributions. 
Further the embedding of the detections due to clusters into 
the LSS is taken into account in an extended model. In Sect. 
[^GEV and order statistics are introduced briefly as well as 
their connection to cosmology. The model predictions are 
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then compared to numerical Wbody simulations in Sect. 
Also the model presented in Maturi et al. (2011) is com¬ 
pared to the simulation. In Sect. we discuss the optimal 
splitting of a wide field survey into sub-surveys in order to 
project out the desired objects, i.e. clusters. The constraints 
on cosmology are finally discussed in Sect. 


2 MODELLING NUMBER COUNTS 


In this section, we provide an analytical recipe to predict 
the number of peaks in weak lensing maps (For a complete 
review on weak lensing basics see [Bartelmann fc Schnei-| 
der|2001 1. The model includes the contribution from obser¬ 
vational noise (such as shot and shape noise), large scale 
structures and non-linear structures such as galaxy clusters. 
The advantage of having a model predicting all contributions 
at once is that it avoids the difficult task of splitting these 
components and keep those arising from galaxy clusters only 
with the risk of having a sample contaminated by spurious 
detections and biased constraints of the mass function. 


2.1 Weak gravitational lensing 

The properties of an isolated lens are given by its lensing 
potential 


0 


( 1 ) 


where $ is the Newtonian potential and Da, Dd and Dda 
are the angular diameter distances between observer and 
source, observer and lens and lens and source respectively. 
The integral is carried out along the line-of-sight, with s 
being the physical distance to the source, ip relates angular 
positions of the source /3 to positions of its images 0 via the 
lens equation (3 — 6 ~ a, where a = Vi/> is the deflection 
angle. For small deflections, the lens equation can be locally 
linearised, introducing a linear mapping described by the 
Jacobian 


V -92 


-92 \ 
l + 9ij ’ 


( 2 ) 


where gt = ji/{1 — k) is the reduced shear and k = ^Aip is 
the convergence. The factor (1 — k) describes the isotropic 
magnification of the image, while the complex shear, g, de¬ 
scribes its distortions. Only the reduced shear g can be ob¬ 
tained by measurements, as the size of the source is un¬ 
known. However, in the weak lensing regime, for which 
«: < 1, 7 « g. 


2.2 Lensing signal 


Weak lensing peak counts can be described as the sum of 
three contributions: 

i) The signal due to real objects, i.e. galaxy clusters. 
Their intrinsic abundance is given by the mass function 
([Press fc Schechter||1974| [Jenkins et al.|2001| [Sheth fc Tor-[ 


men|199'9 1 and the expected density profile. 

ii) The signal due to chance projections of the LSS 
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which is described by the projected two-dimensional power 
spectrum is given by |Limber|1953 1 


P.i£) = 


4cf^ 




a?{w) \fK{w)’ 


(3) 


where mu is the comoving distance to the horizon and W(w) 
is a weight function including the line of sight integral over 
the distribution of background sources G{w), 


fKjw' - w) 


W{w) = I dw'Giw') 


(4) 


In this work the distribution G{w) is taken from Benjamin 
|et al.| p007| ) with parameters (a, 6, c) = (0.748,3.932,0.8). 
Note that this parameter choice implies a distribution with 
mean redshift (z) « 1. As Eq. ([^ gives the total conver¬ 
gence, the tangential component used in Q is Pj^ = Pk/’^ 
due to isotropy. 

iii) The observational noise is given by the intrinsic el- 
lipticity and the finite number of galaxies used to measure 
the shear. This contribution has a white power spectrum: 


P. = 


2ns 


(5) 


with the number density of source galaxies Ug = 40 arcmin“^ 
and ellipticity variance = 0.3 (When comparison with 
the ray-tracing simulation Ug = 30 arcmin“^ and = 0.25 
are used.) 

As the contributions by LSS and the observational noise 
are described by Gaussian random helds, the total contribu¬ 
tion to the signal is given by the sum of their respective 
power spectra. 


P{1) = P-,,{1) + P,. 


2.3 Optimal filtering 


( 6 ) 


A method to measure the signal of weak gravitational lensing 
is the aperture mass ( |Schneid^|1996[ ), which is a weighted 
average of the tangential shear at the position 0: 


A{0)^ f d"0'7t(0',0)g(|0'-0|). 

jRo 


(7) 


Since we are interested in the signal of non-linear structures 
and therefore in the exponential tail of the mass function 
this hlter should maximize the 5'/A^-ratio of non-linear struc¬ 
tures. A linear filter constructed for this purpose was intro¬ 


duced by Maturi et al. (20051 and reads in Fourier space: 


Q(i) = a 


P{£)’ 


where a = 


/ 

Mo 




IfWl 

P{£) 


( 8 ) 


Here P{£) is given by equation § and r is the expected 
shear profile of halos, here we assume NFW halos (Navarro 
et al.||1996[) for wh ich the shear prohle is given analytically 
I Bartelmann|1996 1. The variance of the aperture mass esti¬ 
mate in polar coordinates is 


OO 

/ 

0 




(9) 


When applied to ellipticity catalogues of galaxies the aper¬ 
ture mass ([^ can be approximated by 


Aie) = -Y,^tA0)Q{Wi-s\), 

r) f ^ 


( 10 ) 


where n is the number density of galaxies and et,i is the tan¬ 
gential ellipticity of the i-th galaxy with respect to the posi¬ 
tion 0 which provides an estimate for the tangential shear. 
Q is the filter function in real space given by the Fourier 
transform of Eq. (Isll . 


2.4 Number counts for linear and non-linear 
structures 

The abundance of weak lensing peaks can be described by: 

(1) The contribution by non-linear structures which de¬ 
pends on the mass function and their shear profile. Due to 
shot noise from the discrete positios of background galax¬ 
ies and their intrinsic ellipticities, its contribution to the 
aperture mass A given by Eq. ( |10[ | will scatter around its 
expectation value A{M), where M is the halo mass, with the 
variance given by Eq. Thus halos with mass M will not 
have a unique signal amplitude, but a probability to pro¬ 
duce a certain signal p{A\M), which can be modeled as a 
Gaussian with width a a - It has been shown by |Bartelmann| 


et al. (20021 that the probability for the 5'/N-ratio= S to 


exceed some threshold 5th is given by 


p{S > 5th|-M, z) = -erfc 


5(M, z) — 5th 


V2 


( 11 ) 


Thus, the number of lenses in a redshift interval dz and mass 
interval dM is given by 

d^nni(5 > 5th, M, z) 


= p{S > Sth\M, z) n{M, z), 


( 12 ) 


dzdM 

where n(M, z) is the mass function (in this work we present 
the difference between the [Press fc Schechter|1974| | Jenkins] 
|et al.|2001| and |Sheth &: Tormen|1999 mass functions). Inte¬ 
grating over the volume dV and dM yields the total number 


of weak lensing peaks caused by clusters (see Maturi et al. 

[MT| . 

(2) The contribution by the LSS and the instrumental 
noise which can both be described by a Gaussian random 
field. In fact, the convergence field is not perfectly Gaussian, 
however if the applied filter is chosen to be large enough, 


only small deviations from Gaussianity remain (Maturi et al. 
[2009] ) 

The number density of detections above some signal 
threshold 5 = S/Nua = xua is given by 


niin(> 5) = 


4-\/27r^/^ 


uo 


exp 


Sf_ 

2al 


(13) 


where are the spectral moments generated by the LSS 


and the noise: 


OO 

/ 


a? = / ^"^^P(f)|Q(f)|". 


^ ' 2n 

0 


(14) 


Finally the total number density of detections is given 
as the sum 


ndet(> x) = nni(> x) + niin(> a;). 


(15) 
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In this work we distinguish between negative and positive 
detections. Positive peaks are caused by over-densities and 
thus also by non-linear objects, hence both terms of Eq. (151 
have to be taken into account. In contrast the negative detec¬ 
tions evolve linearly only and are therefore described by the 
LSS term only. Note again that we use the dehnition of de¬ 
tections ( Maturi et al.|2009 l rather than peaks. A detection 
is thus given by a continuous region in a leasing catalogue 
above a given threshold. However, for our application these 
definitions tend to coincide, as we are only interested in the 
highest peaks. 


2.5 Correction term 


The theoretical model for the positive detections according 
to Eq. ( |15[ ) does not reproduce the abundance found in N- 
body simulations (see Figure 5 and Sect.|^, because the two 
contributions entering in Eq. (151 are not independent. In 
fact Eq. (12 I implicitly assumes that the clusters are placed 


on the mean signal (which is zero) of the background caused 
by the LSS. However, the LSS randomly shifts the expected 
signal for a halo around its mean value. We account for this 
effect by thinking of rini as a cumulative distribution, i.e. the 
number of objects above a certain threshold. This random 
shift will move peaks from lower to higher signals and reverse 
with equal probability. However, as there are more objects at 
a lower signal, this shift will effectively increase the number 
of objects at higher signals. To see this, consider some S/N 
threshold (5'/A'^)th and a Gaussian PDF p(x,ao) with zero 
mean and width 


-2 

O-Q 


O-Q 


(16) 


where (Tq is the variance due to the LSS together with the 
noise contribution and cta is the variance of the aperture 
mass due to the noise. At a certain S/N a small fraction 
of detections moves randomly around its expected value to 
larger or smaller values. Therefore the corrected number 
counts are, again with S/N = x, 


nX’^(> x) = nni(> x) -I- gain(a;) — loss( 2 ;), 


(17) 


where loss(a:) is the number of detections which are removed 
from the S/N bin because their amplitude became too large 
or small due to the scatter, while gain(x) is the number of 
detections which are moved from other bins into the consid¬ 
ered bin. Thus, more explicitly we have: 


X 

x) = nni(> x)- j dinni(> x)p{x - x\ do) 

^min 

oo 

-I- J dfnni(> x)p{x — x;ao), 


(18) 


where the prime denotes the derivative with respect to x. 
The additional minus sign is due to the fact that the number 
density at the threshold x is needed, which is given by the 
negative derivative of nni(> x). The symmetry is now broken 
due to the fact that n^i is monotonically decreasing. In the 
upper panel of [Figure I] we compare the corrected model 
(black line, Eq. (13l + (18l) and the uncorrected (blue line, 
Eq. ( |15| )) mode with a raytracing simulation (see sect.j^for 
details). 



S/N 



Figure 1. Top-. Comparison of the total counts above some S/N- 
ratio between model and simulation. The blue dotted curve shows 
the uncorrected model rini-l-niin, while the black solid curve shows 
-b uiin from Eq. | |18[ l. Bottom: The relative impact of the 
correction as a function of S/N. 


The bottom panel of [Figure Ij shows the variation due 
to the correction. The increase is roughly 40% for S/N > 6, 
but negligible for S/N < 3, where spurious detections due 
to LSS and noise contribute strongly to the WL field, thus 
the impact of the correction is small as it does not influence 
the linear regime. Altogether the correction is small at low 
S/N but becomes large in the relevant S/N ranges. 

We finally note remaining discrepancies. Especially for 
S/N £ [3,5], which is the regime where the contribution of 
non-linear and linear counts are comparable. In this regime 
the detections caused by the linear contribution and the non¬ 
linear structures will start to blend affecting the statistics in 
a way which cannot be modeled analytically. A possible way 
out would be the creation of a mock WL map based on a 
Gaussian random fields with power spectrum P{k) given by 
Eq. (§ and adding halos according to their mass function 
and signal prohle. A similar approach has been followed by 


Lin & Kilbinger (2014 \. However, this regime is not of inter¬ 


est for our study as we consider only the peaks with extreme 
values (with large positive S/N) far away from this regime. 
Additionally we are interested in an analytic prescription to 
evaluate the extreme value statistics. 

Furthermore we neglected the influence of halo triaxial- 
ity (see |Oguri fc Blandford||2009 and Hamana et al.|[2012 |. 
This effect will introduce an additional scatter into the de¬ 
tection amplitude. 
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Figure 2. The CDF according to Eq. ( |19[ l of the highest positive 
detections assuming a survey area A = 1 deg^ for different mass 
functions, namely the Sheth-Tormen (solid black line), Jenkins 
(dashed red line) and the Press-Schechter mass functions (dashed 
blue line). 



Figure 3. The CDFs of the first 22 orders, i.e. F(n— 21 ) to ^(JV)i 
of the positive detections for a survey area of A = 1 deg^. Also 
shown is the corresponding GEV distribution as a dashed line. 


3 GENERAL EXTREME VALUE AND ORDER 
STATISTICS 


This section will review the basic concepts of general ex¬ 
treme value (GEV) statistics and order statistics and how 
they can be calculated from the fiducial distribution given 
by the weak lensing counts given by Eq. (151 and (18 1 . 


3.1 General extreme value statistics 


The statistics of the maxima and minima of independent and 
identically distribnted random variables is called extreme 
value statistics | |GumbeI||1958[ ). Consider independent and 
identically distributed random variables Xi and their block 


maximum = max(Vi,..., V„). As shown by Jenkinson 
( |1955[ ) and others, the limiting cumulative distribution func¬ 
tion (CDE) of M„ = X for n —>■ oo is given by 


^-y,a,p{x) — 


exp(^- [l + 7 (^)]' 
exp (exp 


for 7 7 ^ 0 
for 7 = 0 


(19) 


where 7 is the shape parameter, while a and /? are the so 
called location parameters. 7 = 0 corresponds to a Gumbel 
distribution, 7 < 0 to a Weibull distribution and 7 > 0 to a 
Frechet distribution. The corresponding probability density 
function (PDF) can be obtained by differentiating the CDF 
and is denoted with p.y,a,; 3 - To evaluate the GEV parameters 
we consider the weak lensing signal S/N = a; as a random 
variable. The CDF of having the largest signal is the prob¬ 
ability of finding no weak lensing peak larger than x. We 
denote this probability Po{x). Neglecting correlations, i.e. 
the clustering of clusters, (Davis et al. 20111, Poix) is de¬ 


scribed by a Poisson distribution for zero occurrence. Thus 


X 

'i= prob(Xmax X^ = J' (^max)dXn 


= Po{x) = 


exp(—A) 


( 20 ) 


k\ 


exp(-A), 


where A is the expected number of peaks in some patch of the 
sky having a signal larger than x. Therefore A = An^etix), 
where A is the snrvey area and nuet (i) is the nnmber density 


of detections given by (151 or (181 for the old and new model 
respectively. Now the parameters a, /?, 7 from Eq. (191 can 


be derived for the specific peak number density ndet(a:) by 
applying a Taylor expansion around the maximum of the 
corresponding PDFs and equating coefficients. This yields 


7 = n{xo)A-l, P = - 


(1 + 7 )^+^ 

n'{xo)A 


a = xo —— ( 1 + 7 ) ^ —1, 
7 

( 21 ) 


where the prime denotes the derivative with respect to x and 
xo is the most likely signal which can be found by running 
a root search on 


n"{xo) - A{n'{xo))^ = 0, 


( 22 ) 


as a result of the constraint Pq (xq) = 0 (see appendix A.l 


for a derivation). The distribution of the largest weak lensing 
signals can thus be described by the three GEV parameters 
a, P, 7 which can be directly obtained from the number 
counts. 

As an example we show in |Figure 2| the GDFs resulting 
from Eq. (20 1 for a one square degree survey in the fiducial 


cosmology, introduc ed in Sect. [Tj a nd mass fu nctions: [Press 
& Schechter (19741, Jenkins et al. (20011 and Sheth & Tor- 


(19991 in blue, red and black, respectively. It can be 


seen that the choice of the mass function has only a small 
impact on the distribution of the highest lensing signals ex¬ 
cept for the PS mass function. 


3.2 Order statistics 


As described in Waizmann et al. (20131 we now introduce 
the basic concepts of order statistics, describing the statistics 
not of the maxima or minima of a distribution, but of their 
i-th highest or lowest values. 

Suppose Xi,X 2 ,...,X„ is a sample of random num¬ 
bers drawn from some continuous population with a PDF 
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Figure 4. Joint PDFs for different combinations of rank orders of the positive peaks. The calculations are carried out with a survey 
area of 1 deg^. The color ranges from the minimum (blue) to the maximum (red) of the distribution. 


f(x) and a corresponding CDF F{x). The random num¬ 
bers can now be ordered by magnitude and re-labelled 
V( 1 ),X( 2 ), such that V(i) sj Now, Wi be¬ 

longing to X, i.e. the event x < X(i) < a; -|- iJo: we have a 
number of (i — 1) such that sC x. Here exactly one 
Xk lies in the mentioned interval and the remaining (n — i) 
Xk have Xk > a; -|- Sx. Since n observations (random num¬ 
bers) belong to the set, they can be arranged in several ways. 
The total number of such arrangements is 

Tl\ 

A{n,i)= (23) 


with a corresponding probability of having a certain realiza¬ 
tion 


prob(a;) = [T’(a;)]' ^[^(x-I-(5a;) — T’(a:)][l — T’(x)]" (24) 


Here 5x denotes the width of the interval containing only 
one element. For large n, Sx is small and thus terms 0(Sx^} 
can be neglected yielding 




prob (x < X(i) ^ X -|- Sx) 

lim -^ 

Sx^O Sx 


= A{n,i)[Fix)r-^[l-F{x)r-^f{x). 


(25) 


Integrating this equation yields the corresponding CDF, 
which is given by 

= E Q ( 26 ) 


Thus the distribution of the largest, i.e. the extreme value, 
or smallest value is 


F(n)ix) = [F’(x)]" and ^( 1 ) = 1 - [1 - F{x)Y 


(27) 


respectively. This means that, given a set of random vari¬ 
ables Xi, i € [l,n-] drawn from a distribution F{x), the 
smallest Xi can be described as a random number drawn 
from while the largest value is a random variable from 
F’(„)(x). As discussed in Sect. 3.1 for large n both cases can 


be described by Eq. (191. 


This procedure can be extended to joint distributions 
of several orders. Given two ordered random variables X^^), 
such that 1 ^ r < s ^ n and x < y their joint PDF is 


f(r)(s) {x, y) 


(r — l)!(s — r — l)!(n — s)! 

X [F{x)Y-^[F{y) - Fix)r-^-\l - E(y)]"-^ 

(28) 

while the joint CDF is given by 

n 3 I 

F(r)(s){x,y) EE j!(j _*)!(„_ j)! ^gg) 

X [F{x)Y[F{y)-F{x)Y-Yl-F{y)Y^-F 


[Figure 3| shows the Erst 22 orders in a one square degree field 
for the positive detections. The expected number of obser¬ 
vations is therefore N = ndet(xmin)A, with the detection 
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threshold Xmin- Note that we need to introduce the cut-off 
Xmin because ridet cannot be converted into a distribution 


function as it is not monotonically decreasing (Maturi et al. 
20091. However, by setting a boundary, x-min, we restrict our¬ 


selves to the region where the distribution is monotonically 
decreasing. The actual choice of a;inin is not crucial, as long 
as it converts ridet (*) into a monotonically decreasing func¬ 
tion for all X > Xmin and as long as it is small with respect 
to the expected values x given by the resulting PDF. We 
thus write the CDF of weak lensing peak counts as 


F{x) = 1 - 


ridet (x) 
ridet (Xmin 


(30) 


which enters in Eq. (291 and (261. As a first result it can be 


seen that the CDF steepens with increasing order and thus 
the constraining power of a single measurement increases. 
This is due to the fact that for detections with lower S/N 
more objects occur in the corresponding signal bin. We will 
call the detection with he highest signal the first order, the 
detection with the second highest signal the second order 
and so on. 

[Figure 4| shows the joint probability distribution of dif¬ 
ferent combinations of orders, as denoted in the individual 
panels. For example we compare the first with the third or¬ 
der in the upper left plot of [Figure 4| which quantify the 
probability of finding the largest weak lensing signal at a 
certain S/N, while finding the third largest signal at an¬ 
other S/N. All calculations were done with A — I deg® and 
Xmin = 1.5, where Xmin is again the threshold introduced 
before. One can see that the closer the orders are to each 
other the higher their correlation is. In addition, higher or¬ 
ders confine the PDF to a substantially smaller S/N region. 
This reflects the fact that the PDF steepens for higher order 
combinations, restricting the PDFs to a smaller range in the 


S/N-pla,ne as already shown in Figure 3 


4 CONFRONTING THEORETICAL 

PREDICTIONS WITH SIMULATIONS 

4.1 A-body numerical simulation 

We test our prediction outlined in Section and by pro¬ 
cessing a mock ellipticity catalogue derived by a ray-tracing 
A-body numerical simulation carried out with GADGET-2 
I Springel|2005 I and presented by Borgani et al. (20041. The 
simulation represents a standard ACDM model including 
dark matter as well as additional baryonic physics for the 
gas component modeled by radiative cooling, star formation 
and supernova feedback in the form of galactic winds. The 
base cosmological parameters are Hm = 0.3, Ha = 0.7 and 
Hi, = 0.04. The Hubble constant is Ho = lOO/ikms”^ Mpc“^ 
with h — 0.7 and the linear power spectrum of the matter- 
density fluctuations is normalised to as = 0.8. The simu¬ 
lated box has a side length of 192/i“^Mpc, containing 480® 
dark-matter particles with a mass of 6.6 x 10^h~^M q each 
and an initially equal number of gas particles with a mass 
of 8.9 X 10^h~^M q each. 

As discussed in [Pace et al.| (|2007|) and [Maturi et al.| 


120091 backward light cones have been constructed using the 


simulation by taking snapshots at different redshifts ranging 
from z = 1 to 2 = 0. In order to avoid repetitions of the 
same structures in snapshots at different redshifts, these are 


randomly shifted and rotated. The light cone has been sliced 
into thick planes whose particles are subsequently projected 
with a triangular-cloud scheme on lens planes perpendicular 
to the line-of-sight. Finally a bundle of 2048 x 2048 light rays 
is traced through the light cone starting from the observer 
into directions on a grid of 4.9 degrees on each side. From 
this ray tracing simulation effective convergence and shear 
maps are obtained. These are used to lens a background 
source population of galaxies at 2 = 1 with a number density 
of rig = 30 arcmin“® and random intrinsic ellipticities drawn 
from 


exp[(l - e®)/o-® 
7rcr2[exp(l/o-|) - 1] ’ 


(31) 


with (Je = 0.25. The mock catalog of galaxy ellipticities is 
then analysed with optimal filter described in Sect. [2.3] on a 
grid of 512 x 512 positions covering the entire field-of-view. 
To account for the survey geometry, the pixelized shear map 
and the discretization of the sky, a filter functions W needs 
to be introduced in the spectral moments (Eq.(14l when 
evaluating the number counts model from the simulation. 
The result of the filtering procedure is 


^2 = j ^"^^P(fc)W®|Q(fc)|®, (32) 

where Q is the optimal filter introduced in Eq. ([^ and the 
filter function VF is a product of the following terms: a high- 
pass filter, suppressing scales larger than the side length of 
the field of view. 


Wt®(A:) = exp 


(33) 


a low-pass filter, filtering out signatures on scales smaller 
than the average separation of sources. 


Wg{k) = exp 



(34) 


and another low-pass filter including the finite size of the 
pixels 




2y/TT j f kd 

pix \ 

fcdpix ^ V \/^ / ’ 


(35) 


where kf = 27r/Lf, fcg = 27ry^fi^,dpix = 4.9°/512 and Ji 
is a cylindrical Bessel function of order one, representing a 
circular step function of the size of a pixel. 


4.2 Estimation of the CDF 


The detections are defined by the up-crossing criteria dis¬ 
cussed in Maturi et al. ( 2009| ), i.e. the resulting weak lensing 
maps are sliced into different S/A-ratios. Then the number 
of detections is counted above every threshold. A detection 
is defined as a continuous group of pixels above the given 
threshold, which are linked if they are closer to each other 
than a certain linking length. In this case the latter is 1.42 
pixels, so that also pixels diagonal with respect to each other 
are linked together. 

In order to study the order statistics of the WL maps 
in the numerical simulation, the field has to be divided into 
A smaller sub-fields, each with area Asub. The sample of 
the i-th highest value drawn from each sub-field represents 
the order statistic of order i corresponding to a survey area 
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Figure 5. Comparing out analytic predictions with a ray-tracing A^-body si mulat ion (grey crosses with errorbars). The survey area, A, 
is 4.9 X 4.9/24 deg^ and the cosmological base parameters discussed in Sect. [4T] are used. Top: Prediction for the negative detections 
(black line), i.e. linear structures only (Eq. Bottom: The model discussed in |Maturi et al.| l |2011[ | (see Eq. ( |15| |) is shown in blue 

(dashed), while the new model including the correction from Eq. ( |18[ l is shown in black. 


Asuh’ The CDF of the z-th order from the simulation can be 
estimated by 

^w(4)) = ^J’ (36) 

where is the ordered (from small to large values) vector 
of the i-th highest value in each subfield Aaub- As all are 
drawn from the same distribution, /(i)(x), their respective 
error is given by the width of the distribution. In Sect. |4.3| 
and |4.4| we compare the order statistics of both, positive and 
negative S/N, with the simulation. By splitting the WL map 
into 24 tiles. Note that the error bars are rather confidence 
regions which are given by the width A(i) of the PDF of the 
i-th order statistics. We define this to be 

A{i) = Qso — < 320 , (37) 

where Qp is the quantile of the distribution with p = 100 x 
CDF(Qp), i.e. the inverse of CDF. The value x^^j of the 
simulation is than given by x|^j = x^^j ± A(i)/2. 


4.3 Negative detections 

The first row of [Figure 5| shows the model prediction for 
the first three orders together with the distribution result¬ 
ing from the numerical simulation calculated from Eq. ( |36[ ). 
Clearly the prediction lies well within the error bars. The 


steepening and shifting to smaller values of the order statis¬ 
tic in the simulation are in good agreement with the theo¬ 
retical prediction. We stress the fact that [Figure 5| shows a 
prediction and not a fit of the model to the simulated data. 


4.4 Positive detections 


In the second row of [Figure^ we repeat the same exercise for 
the positive detections. Again, the model is in good agree¬ 
ment with the simulation. Note that for the first order the 
CDF from the simulation cannot reach unity because there 
are two peaks at very high S/N, namely S/N = 13 and 
S/N = 17, which are not shown for easier comparison of 
the different orders. In this Figure two models are shown: 
The dashed blue line corresponds to the order statistics cal¬ 
culated from the uncorrected version of the number counts, 
while the black line includes the correction acting on the 
detections caused by non-linear structures. This clearly il¬ 
lustrates that the contribution of the LSS is relevant also 
for the detections with the highest S/N ratios. Ignoring the 
contribution would lead to biased cosmological parameters, 
even though the detections at large S/N are only due to 
non-linear structures. We further notice that the difference 
induced by the correction is much larger than the one due 


to the use of different mass functions (Figure 2|. 
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5 OPTIMIZING THE AREA OF THE 
SUB-FIELDS 


The previous analysis has been performed by splitting the 
simulated field (4.9 x 4.9 deg^) into N sub-fields. So far 
the choice was arbitrary and we used a size of 1 deg^ for 
simplicity. Thus, the question is whether this can be done 
in a more objective and optimal way to do it. 

Since the goal is to estimate some distribution on a large 
number of sub-fields is preferred as the sampling improves 
and thus the overall constraining power on cosmological pa¬ 
rameters. However, increasing the number of sub-fields au¬ 
tomatically decreases the size of each sub-field given a fixed 
total survey area, A. Therefore the strongest lenses in each 
sub-field will have smaller signals, diluting the presence of 
those with larger amplitude and mostly related to non-linear 
structures. Furthermore the sensitivity on cosmological pa¬ 
rameters shrinks due to the dilution of non-linear structures 
and LSS due to the noise. The probability of finding a signal 
above some threshold tsn is 

OO 

Prob[U„;A] = J pNix)dx, (38) 

tan 


where PN(a:) denotes the PDF of the order statistics corre¬ 
sponding to a sub-field size of Asub = A/N. Again, increas¬ 
ing N will lead to a better sampling of the distribution, but 
it will also lead to a lower probability of finding a detection 
due to a cluster in each sub-field. The expected number of 
objects above this threshold is thus 

(N) = NPToh[tsn;N]. (39) 


The threshold chosen such that the largest signals are 
likely to be due to clusters. Eq. (381 is the probability that 
a halo with S/N > tan can be found in a certain area, con¬ 
versely, 


q = I — Prob [tan, N] = 1—p (40) 

is the probability of not finding a halo. However, aiming at 
galaxy clusters the largest signal in each sub-field is required 
to be due to the non-linear contribution. Therefore no object 
should be below the threshold tan which marks a boundary 
between non-linear and linear counts, thus 


qN ^ 1, 


(41) 


such that in each sub-field the corresponding largest signal 


occurs above the threshold. Eq. (411 sets a maximum for N. 


Combining this with the goal to sample the distribution as 
well as possible the maximum N can be used as the number 
of sub-fields. As a result we present in |Figure 6| the number 
of sub-fields as a function of the threshold tan and the total 
survey area. As expected, increasing the threshold lowers the 
maximum number N, since there are fewer objects above 
higher thresholds. Accordingly with increasing total survey 
area the N also increases as more data are available. As 
the total survey area is fixed by the experiment the only 
free parameter left is tan, which should be chosen such that 
everything above the threshold will most likely be a halo 
and not due to LSS. A reasonable value for the threshold is 
found to be tan « 5 — 6, which is the value at which well 
observable galaxy clusters show themselves. 



survey area 


Figure 6. Dependence on the optimal number of sub-fields on 
the threshold tan and the total survey area. 



Figure 7. The relative sample bias as a function of the number 
of sub-fields for a fixed sub-field size A. and trg were fitted 
separately. 


6 CONSTRAINING COSMOLOGY 


In this section we use the method outlined in Sect. |2] and 
[^to constrain cosmological parameters. Before applying the 
method to mock data we briefly discuss the impact on indi¬ 
vidual cosmological parameters on the extreme value statis¬ 
tics. Afterwards we carry out a likelihood analysis using a 
Euclid-like survey with A — 15000 deg^. We mimic a source 
redshift distribution with a mean redshift of approximately 
unity, using the parametrisation given in [Benjamin et ah] 
(20071 {a,b,c) = (0.748,3.932,0.8) and a source number 


density rig = 40arcmin“ . We investigate two scenarios, 
one completely ignoring redshift information and one divid¬ 
ing the survey into two redsift bins. In the first case we 
use N = 4574 sub-fields and in the second case the follow¬ 
ing number of sub-fields: Wow = 2858 and Ahigh = 2059 
according to the criterion given in Sect. Mock data are 
created by sampling N random numbers drawn from the dis¬ 
tribution of the largest signal given by Eq. \22\ . Note that 
small survey sizes, allowing a sampling of A ~ 500 sub¬ 
fields, can create a sample bias when deriving the cosmo¬ 
logical parameters because of under-sampling of the fiducial 
distribution and its non-Gaussianity. This effect can be re- 
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Figure 8. Relative deviations A in percent of quantiles of the extreme value distribution (1st order) as a function of cosmological 
parameters. Quantiles used are 5% (solid blue) 50% (dashed black) and 95% (dotted red). The fiducial model is given by the base 
parameters described in [section 1| Top -, positive detections. Bottom: negative detections. 


duced by using the average likelihood, i.e. by averaging over 
many realisations of Eq. |44| This will increase the uncer¬ 
tainties in the derived parameters. In [Figure 7| the relative 
bias {9fit — 9fLd)/9fid, with 9 being the parameter is shown 
as a function of the number of sub-fields N. It can be seen 
that the bias increases at low N. At higher N it is close to 
zero and approximately constant. The systematic offset just 
arises from the finite resolution of the likelihood and can be 
in principle reduced to zero by introducing more grid points 
in parameter space. Note that the plotted parameters were 
not fitted simultaneously. However, for our application the 
number of sub-fields is much larger and this sample bias will 
not pose any problem. 


6.1 Impact on cosmology 

We take a first look at the dependence on cosmology of the 
order statistic ignoring any redshift information as a starting 
point. 

[Figure ^ shows the dependence of quantiles Qp, i.e. the 
inverse of the CDF at a given value p of the extreme value 
distribution (1st order) of both positive and negative peaks 
as a function of cosmological parameters (for example Qso 
corresponds to p = 0.5 and therefore 50% probability occurs 
below Qso). We show the relative deviation from the fiducial 


model 

A = {Qp- gf )/0“, (42) 

and for which we assumed a fiat universe. 

The sensitivities to Dm and as are comparable and very 
strong in case of the positive detections. This also shows the 
degeneracy between the two parameters which is typical for 
estimates based on weak lensing. The shape of the distri¬ 
bution changes slightly when varying the parameters as the 
curves of different quantiles have different slopes. This effect 
is due to the highly non-linear dependence of the positive 
detections on those parameters. This behaviour cannot be 
seen in the bottom row, where the negative detections are 
shown. Because of the Gaussian nature of the peaks in this 
case, the overall shape of the extreme value distribution does 
not vary significantly when changing cosmological parame¬ 
ters. We also note that the sensitivity is much stronger in 
case of the positive detections, causing changes up to 12% 
in case of flm and ag in the considered region of parameter 
space, while the negative detections only give rise to changes 
up to 4%. 

The dependence on the equation of state parameter is 
much weaker than the dependence on parameters related 
to structure formation. Volume effects are less important, 
compared to the dependence on Dm and ag. The increase of 
the relative difference for smaller wq (less negative) is mainly 
due to the modifications of the exponential tail of the mass 
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Figure 9. Influence of the source number number density Ug on 
the CDF for fixed area A. 


function. It can again be seen that the negative detections 
show a weak sensitivity compared to the positive peaks with 
respect to wq. This agrees with the fact that volume effects 
are not the dominating factor. 

Note that when adopting a time dependent equation 
of state w{a) = wo + Wa.{l — a), the quantiles nearly stay 
constant in the parameter range. In Waizmann et al. (20131 
and Pace et al. (20101 it has already been shown that the 
most massive objects are only affected at higher redshifts 
when adopting a simple linear model for w{a). For halos 
at high redshift the lensing signal shrinks and thus those 
objects are not captured. 

Additionally, the number of background galaxies rig 
shrinks when going to higher redshifts, thus the S/N-ratio 
drops as the shot noise increases. Thus, even though the 
abundance of very massive halos will change (see [Waizmann 


et al.|[ 201 ^, the most extreme weak lenses are not affected 


that much as those massive halos appear at higher redshift, 
where the S/N-ratio is lower in the first place. Therefore the 
dependence on rua is very weak as its impact becomes impor¬ 
tant only at higher redshift when keeping wo constant. Note 
that the impact of cosmological parameters also depends 
on the adopted source redshift distribution, for example a 
higher sensitivity in wo can be reached if the sources popu¬ 
late higher redshifts. 

We finally note that the corresponding GEV parame¬ 
ters Q, /3 and 7 depend nearly linearly on the cosmological 
parameters except for wo, allowing the use of simple schemes 
to quickly evaluate the likelihood. 


6.2 The likelihood 


For calculating the likelihood a mock data set is created 
by sampling N random numbers from the fiducial model. 


i.e. from the derivative of Eq. (20 1 Pa,p,-y{x) evaluated at 
the fiducial cosmology. The components of the data vector 
D are denoted {D)i = A;.Those N numbers are the highest 
5'/Al-ratios in the N sub-fields with area A = Atot/N, i.e. we 
fix Atot = 15000 deg^ and calculate the optimal number of 


sub-surveys N via the condition (411. Having fixed Atot and 


N, the mock data is sampled from the resulting distribution 
( | 20 [ ) for both the positive and negative peaks. 

The likelihood is now the joint conditional probability 


Figure 10. Likelihood for a Euclid like survey with A = 15000 
deg^. and without redshift information. We show confidence re¬ 
gions corresponding to la, 2a and 3a. Positive and negative 
counts are used, i.e. j = ± in Eq. | |44| l. The fit was carried out 
with ag and only. The fiducial cosmology is marked with a 
dot. 


for having a set of parameters, collected in the parameter 
vector 0 , given a data set D: 

N 

£(©lZ?) = nPc,/3,^,>(©|Ai), (43) 

where we already assumed a flat prior p{&) and the evidence 
p{D) was already included into the normalization. The log- 
likelihood L is 

N 

^(©) = E E log , (44) 

i = l J = ± 

where the plus and minus signs denote positive and negative 
peaks respectively. Again flatness is assumed, i.e. Ha = 1 — 
Hm, while parameters like h and Hb are fixed in the first 
place, as their influence is negligible, thus yielding only an 
overall factor. Since the individual terms Pa, 0 ,~i,i,j of the 
likelihood are not Gaussian or have a simple exponential 
form, also the likelihood will of course not assume this simple 
form. Therefore a Fisher approach is not suitable and the 
likelihood has to be sampled directly or via Monte-Garlo 
methods. 

In [Figure 9| we show the influence of the source num¬ 
ber density rig. Clearly, fewer background galaxies lead to a 
lower weak lensing S/N-xatio, which is expected. However, 
the sensitivity of the highest peaks to cosmological parame¬ 
ters is due to the mass function. Thus the shape of the like¬ 
lihood will be conserved. Nonetheless, as the S/N decreases 
if rig is decreased, the expected S/N of the highest peak will 
also decrease if the sub-field area is kept fixed. For example, 
if we are given two surveys with the same total survey area 
Atot but different source densities rig, the possible number of 
sub-fields N will also differ. Especially if rigi > ng 2 we have 
Al > A 2 . As the sub-fields are considered to be indepen¬ 
dent, the likelihood scales with (;J^e to the Poissonian 

nature of the statistical error. In the survey discussed in 
this work the uncertainties of the cosmological parameters 
for rig = 40 deg“^ will pick up a factor •y/A(40)/A(30) « 1.2 
if we calculate the likelihood with rig = 30deg“^. 

When including the two redshift bins, their individual 
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contributions are included in Eq. (44 1 as an additional sum 


over all redshift bins. Furthermore the upper limit of the 
summation N is bin-dependent. The likelihood therefore 
reads Lz{&) = 5I^Lbin(©), where Z/bin(®) is given by Eq. 
(441. 


6.3 Confidence regions 


We show in |Figure 10| the resulting likelihood without using 
redshift information in the flm — os plane. In this case we 
only fit flm and as and keep wo fixed, because its variation 
around the hducial value only marginally affects the results. 
The fiducial cosmology is marked with a point. It can be seen 
that the expected degeneracy between flm and as exists. 
Note that the likelihood is really non-Gaussian even though 
the contours look like ellipses, with contours correspond¬ 
ing to 68.3%, 95.4% and 99.7% confidence. The error con¬ 
tours are nearly orthogonal to those obtained from the CMB 
I Planck-Collaboration et al.|2013 | and thus a joint analysis 
can yield very tight constraints on both parameters. The 
relative uncertainties of both parameters are ~ 10“^. Note 
that, assuming a source background density Ug = 30deg“^ 
the relative uncertainty will pick up a factor of 1.2 as men¬ 
tioned before. Especially for ag this result is competitive. 


for example Wang et al. (20101 showed that a relative un¬ 
certainty of ~ 7 X lO”'” can be achieved by using Baryon 
Acoustic Oscillation measurements. 

So far we neglected any redshift information regard¬ 
ing the detections discarding all the information about the 
time evolution of the universe. Here, we introduce only two 
redshift bins for simplicity in order to trace the evolution 
of structure formation. By introducing redshift boundaries 


when integrating over Eq. (121. We would like to stress that 


here we assume the minimum amount of information when 
including redshift in the analysis, i.e. we do not assume to 
have the redshift of each individual detection, but we sim¬ 
ply perform the peak counts on two shear maps, one for the 
higher and one for the lower redshift bin. This is because a 
redshift estimate can only be associated to those detections 
which have been identihed to be galaxy clusters ( |Bellagamba| 
et al.||2011 Bellagamba et aL]|2012 (. The estimate could be 
given by taking advantage of optical filters or other external 
data. Thus the procedure described in Sect. has to be ap¬ 
plied to each of the two shear maps separately, i.e. by fixing 
the total survey area one can calculate the optimal number 
of sub-fields for each shear map. 

In order to get comparable statistics in the redshift bins, 
the expected number of objects should also be comparable. 
For our purpose we use two redshift bins such that the ex¬ 
pected number of peaks above the threshold tan is roughly 
equal, resulting into redshift bins with 0 < 2 : < 0.35 (named 
’low’) and 0.35 < z (named ’high’). 

We fit Hni, erg and wo simultaneously while setting 
uia = 0. In [Figure ll| we show the resulting likelihoods. On 
the diagonal the maximized likelihoods can be seen, i.e. the 
one dimensional sub-spaces with the other parameters set 
to their best fit value. The lower triangle shows the maxi¬ 
mized constraints on all three parameter permutations. One 
can see that the constraining power in the — ag plane is 
very similar to the one obtained from [Figure fol However, 
Wo can now be constrained allowing an uncertainty below 
0.1 at Icr and thus a relative uncertainty < 10%. We also 


note that the degeneracy between wo and each of the other 
two parameters is very small. 


7 CONCLUSIONS 

The results presented show that applying order statistics 
to the distribution of cosmic shear peaks is a powerful tool 
to constrain cosmological parameters associated with struc¬ 
ture growth and the geometrical evolution of the universe. 
We evaluated the extreme value and order statistics for 
both over-densities and under-densities of weak leasing shear 
maps, we optimized the criterion to split a wide field survey 
into sub-fields in order to optimally sample the distribution 
of extreme values, and we forecasted the constraints on flm, 
ag and wo achievable with wide fields surveys. In the pro¬ 
cess we also improved the analytic recipe provided in [Matun| 
et al. (20091 and Maturi et al.j ( |2011| ) ay accounting for the 
impact of LSS on the larges S/N peaks which cannot be 
neglected even though these are caused by highly non-linear 
structures. A good agreement was found when comparing 
our model with a ray-tracing N-body numerical simulation. 

This has been done with mock data under the assump¬ 
tion that (1) the contribution of noise and LSS in weak leas¬ 
ing shear maps produced with relatively broad kernels are 
well represented by a Gaussian random field and (2) by ig¬ 
noring blending between the peaks mainly caused by noise 
and LSS and those caused by non-linear structures. For a 
Euclid like survey with and without including photometric 
redshifts of the galaxies in our analysis we obtained the fol¬ 
lowing: 


(i) When modelling the statistics of weak leasing peak 
maps the embedding of the contribution from clusters into 
the LSS has to be taken into account. While this might 
appear as counter-intuitive at first glance, but both LSS and 
noise fluctuations (even if occuring at low S/N ratios) are 
very important also in the high signal-to-noise regime. Their 
impact is not negligible, especially when constraining the 
extreme value statistics, and cannot be ignored. The impact 
of the LSS and noise contributions to large S/N ratio peaks 
is of order 15%. For this reason great attention has to be 
given to a detailed description of the data noise which, even 
if of low amplitude, may result in large biases if its statistical 
properties are not well understood. 

(ii) Using extreme value statistics we give an analytic pre¬ 
diction for the largest weak lensing signals expected in wide 
held surveys. These values can also be used to verify if the 
existence of extremely massive clusters such as ‘El gordo’ are 
falsifying ACDM or not. For example, this cluster, which has 
been widely claimed to be ’troublesome’ for ACDM, is not in 
contrast expectations from ACDM ( Waizmann et al.|2012 |. 

(iii) The extreme value statistics applied to peak counts 
has the advantage that the resulting likelihood from which 
cosmological parameters are extracted does not need any 
assumption of the shape of the distribution as it is given 
by the model itself. Furthermore the identiheation of the 
highest peaks is a straightforward task. 

(iv) We evaluated the constraints of the cosmological pa¬ 
rameters achievable with the extreme value statistics applied 
to weak lensing maps. Booth and ag can be well con¬ 
strained below the percent level, additionally the parameter 
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Figure 11. Confidence regions for a survey with A = 15000 deg^) and two redshift bins. Confidence regions shown in the lower triangle 
are again Icr, 2a and 3cr. The plots on the diagonal show one dimensional cuts through the likelihood, while plots in the lower triangle 
show two dimensional cuts and the respective confidence regions. The fiducial cosmology is marked with a dot. 


degeneracy is nearly orthogonal to that obtained with CMB 
measurements. 

(v) The dark energy equation of state can also be con¬ 
strained by splitting the ellipticity catalog in two redshift 
bins so that the absolute error of wo is A(tno) < O-l- This 
result is competitive with other cosmological probes, e.g. 
the constraints obtained by the |Planck-Collaboration et al.| 

Investigating the highest shear peaks in wide weak leas¬ 
ing surveys, the presented method directly reflects the statis¬ 
tics of the most massive objects, i.e. galaxy clusters, avoiding 
any reference to their mass and the danger of using poten¬ 
tially biased mass proxies and scaling relations. Even though 
few assumptions have been used, this method has the po¬ 
tential to provide strong constraints on the underlying cos¬ 
mological model. 
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yields immediately 
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—n'{xo)A = 


P 


1 + 7 


xo — a 


1 -1/7-1 


(A4) 


Due to the maximum constraint the third term of the ex¬ 
pansions has to vanish: 


— 0 — Pq{xq) 

yielding 


1 + 7 = 


1 + 7 


Xo — a 


P 


■ -1/7 


and Eq. (191. Plugging (A 6 | into (A2| yields 
7 = n{xo)A — 1. 


After inserting (A 61 into (A4 1 one can conclude 


P = - 


(i+7r+\ 

n'{xo)A 


(A5) 


(A 6 ) 


(A7) 


(A 8 ) 


Using again (A4| and inserting P from the latter equation 

(A9) 


one ends up with 

a = *0 - ^ [(1 + 7 )”^ - l ] • 


This paper has been typeset from a T 15 X/ DTTiX file prepared 
by the author. 


APPENDIX A: DERIVATION OF THE GEV 
PARAMETERS 


Equating coefficients of the Taylor expansion around the 
peaks of both distributions one finds via 


Po(^o) — G^^a^^{xo') 
the equation 


-n{xo)A = 


1 + 7 


Xo — a 


The second term 

^ 0 (^ 0 ) f^7,a,/3(^o) 
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